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Complementary diffuse and inelastic synchrotron X-ray scattering measurements of lead zirconate- 
titanate single crystals with composition near the morphotropic phase boundary (x=0.475) are 
reported. In the temperature range 293 K < T < 400 K a highly anisotropic quasielastic diffuse 
scattering is observed. Above 400 K this scattering disappears. Its main features can be reproduced 
by model of inhomogeneous lattice deformations caused by inclusions of a tetragonal phase into a 
rhombohedral or monoclinic phase. This observation supports the idea that PZT at its morphotropic 
phase boundary is essentially structurally inhomogeneous. 

PACS numbers: 77.80.Jk, 61.72.Dd 



Lead zirconate-titanate (PbZri-a^Tia^Os, PZT) is one 
of the most technologically important ferroelectrics Q . 
Being widely employed in practice PZT is also a model 
system representing ferroelectric solid solutions with a 
morphotropic phase boundary (MPB). Understanding of 
the mechanisms leading to very high dielectric and piezo- 
electric responses near the MPB is essential for strategic 
design of new and improved materials, particularly the 
ecologically friendly lead- free PZT counterparts Nu- 
merous theoretical and experimental studies performed 
on different lead-containing MPB ferroelectrics produce a 
highly complex and often controversial picture. Many of 
these aspects are covered in a 2006 review by Noheda and 
Cox and references therein. More recent trends are 
briefly highlighted in a 2009 editorial review by Kreisel 
et al. [H. Recently the morphotropic PZT single crys- 
tals of various compositions become available [6], which 
has opened up new experimental possibilities. The first 
single-crystal diffraction experiments [7|, |8[ did not pro- 
vide the final conclusion regarding the true microscopic 
structure of PZT, but allowed to make several important 
conclusions. On the basis of neutron diffraction study of 
morphotropic PZT with x=0.46, it was shown [7] that 
low-temperature monoclinic Cc phase should be ruled 
out as a ground state and a coexistence of rhombohe- 
dral and monoclinic Cm domains should be considered 
instead. High-resolution X-ray diffraction study [H also 
supports a phase coexistence model for that composition. 
By using 2-dimensional single-crystal scattering maps in- 
stead of 1-dimensional powder spectra the authors suc- 
ceeded in resolving otherwise overlapping Bragg reflec- 
tions and demonstrated the presence of more than one 
phase. The idea of phase coexistence was also supported 



by recent studies by anelastic and dielectric spectroscopy 
[oj and neutron powder diffraction [loj. 

A number of times it was also pointed out that mor- 
photropic PZT can be inhomogeneous on the nanoscale. 
This point of view is supported by observation of nano- 
metric contrast fluctuations within micrometer-scale do- 
mains in PZT by transmission electron microscopy 
12]. Twinned nanodomains were also considered as a 
cause for unusual optical isotropy P revealed in the 
tetragonal phase of PZT with x=0.46. From another 
point of view [is] the nanoscale heterogeneity in PZT can 
be connected with regions of short-range correlated mon- 
oclinic ionic displacements which on average produce the 
diffraction pattern compatible with rhombohedral and 
tetragonal symmetry at different sides of the phase dia- 
gram. It was also suggested on the basis of single-crystal 
inelastic X-ray scattering [14] that morphotropic PZT 
has similar to relaxors relaxational-type zone-boundary 
lattice dynamics and shares with them some intrinsic 
nanoscale inhomogeneity. 

A powerful technique of studying structural and other 
types of inhomogeneities in crystals is diffuse scattering 
(DS). Particularly it proves itself useful in the studies of 
structural instabilities in ferroelectrics and related sys- 
tems [HI. To the best of our knowledge, no publica- 
tions have been available in which the diffuse scattering 
was observed or interpreted in PZT except the electron 
diffraction study by Glazer et al jTH . The authors report 
the observation of DS in both Zr-rich and Ti-rich com- 
positions but not in the ones close to the morphotropic 
phase boundary. By using single-crystal synchrotron X- 
ray scattering we show that the diffuse scattering indeed 
exists in morphotropic composition of PZT, evidencing 
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FIG. 1: Diffuse scattering (DS) distributions in the (0 4) zone at room temperature (a) and 683 K (b), DS distribution in the 
(-3 1) zone at room temperature (c), temperature dependence of integrated DS intensity during heating (circles) and coohng 
(triangles) (d). 



its structural heterogeneity. 

Morphotropic PbZri-a^Tia^Oa single crystals with 
PbTiOa content x=0.475 were grown by a top-seeded 
solution method. For the X-ray measurements a stick- 
shaped sample with about 100 x 100 micrometer cross- 
section was prepared by slicing, polishing and subsequent 
etching in HCl acid. Diffraction and diffuse scatter- 
ing measurements were carried out at Swiss-Norwegian 
beamlines at ESRF using KUMA (Oxford Diffraction) 
diffractometer with CCD detector. A locally constructed 
heat blower was used for heating up to 773 K. 

At room temperature we obtained the diffuse scatter- 
ing distributions in the (0 4) and (-3 1) zones (Figs. 
[1^, and [It). They are highly anisotropic and appear to 
resemble the DS shapes in relaxors [16|, [l^- On heat- 
ing this strong DS disappears at temperatures between 
373 K and 423 K and only weak, but also anisotropic 
diffuse halo remains up to 773 K. Characteristic distribu- 
tion of this high-temperature DS is shown in Fig. [TJd. On 
cooling the strong DS reappears in the same temperature 
region, but according to our measurements is systemati- 
cally less intense than before heating. The temperature 
dependence of the DS intensity is presented on Fig. [l]i. 
The intensity points on that plot correspond to the values 
of parameter Iq obtained by data fits described below. 

A log-log plot of diffuse scattering profile along the 
high-intensity direction [0 1-1] in the (0 4) zone is pre- 
sented in Fig. 2. The data for q < 0.02 r.l.u. are spoiled 
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FIG. 2: Diffuse scattering profile along the [0 1 -1] direction 
in the (0 4) Brillouin zone at room temperature together 
with the results of fit by expression I{q) = Ioq~^ . 



by Bragg scattering and not included in the plot. For 
larger q values a power law I{q) = Ioq~^ is observed with 
the power s — 2 ± 0.2. The corresponding fit is shown by 
solid line. For estimation of the temperature dependence 
of DS intensity in Fig. [T]i the value of parameter s was 
rounded up to be 2. 

Diffuse scattering distributions obtained at SNBL in- 
clude both the elastic and inelastic contributions. To 



3 



> 

CD 

E 



rn 500 



100 



> 

CD 

E 



> 

CD 

E 




A 



(1) 



FIG. 3: Inelastic X-ray scattering intensity maps for the di- 
agonal, longitudinal and transverse directions in the (2 0) 
Brillouin zone. 



distinguish the nature of this DS we additionally per- 
formed an inelastic X-ray scattering (IXS) experiment at 
ID28 ESRF beamline. The resolution was about 3 meV. 
The IXS maps along the diagonal (2-h,h,0), longitudinal 
(2+h,0,0) and transverse (2,k,0) directions are presented 
in Fig. [3l These maps apparently demonstrate that the 
maximum signal corresponds to elastic scattering (£^=0) 
for diagonal and longitudinal directions but no elastic line 
is observed for transverse direction. This contrasts with 
the DS in relaxors where temperature-dependent DS ex- 
ists in the transverse direction but is almost absent in 
the longitudinal direction [18]. 

Surprisingly, we do not see any increase of X-ray DS 
intensity near T=663 K, where the cubic-to-tetragonal 
transition takes place. In fact this transition is accom- 
Danied by a strong central peak in Brillouin scattering 
19], most possibly associated with ferroelectric fluctua- 
tions. We do not see such fluctuations by X-rays near 
high-temperature transition and thus do not expect to 
observe them near low-temperature transition. This way 
we see the strong increase of DS intensity below 423 K as 
a sign of developing heterogeneity. A starting point for 
interpreting this heterogeneity can be set up on the basis 
of previous results, that indicate the simultaneous pres- 
ence of multiple phases. We start from an assumption 
that a host phase of specific symmetry contains clusters 
of a different symmetry. In this case the DS can be de- 
scribed by the terms corresponding to the form factor of 
the clusters and the factor describing the impact of these 
clusters on the matrix [2Qi] : 



The intensity is proportional to the number of defect 
centers N^^ Debye- Waller factor e~^^ and a sum of A 
additives. Each additive describes the scattering due to 
particle of orientation a. The shape of particle is rep- 
resented by Fourier transform of corresponding shape 
function Saiq)- A/ represents the difference between the 
structure factors of the host phase and the clusters. The 
term fQA^oc describes the scattering due to elastic defor- 
mations caused in matrix by particles. This latter term 
is assumed to be dominant in comparison with A/ since 
the structures of matrix and particles in PZT are ex- 
pected to be very close. When particles are sufficiently 
small we may also neglect a form factor and focus our 
consideration on elastic deformations. Characteristic 2- 
D distributions of DS are tabulated in Ref. |20| for many 
simple defect symmetries and one may easily find that 
tetragonal defects in cubic matrix will cause DS very 
similar to the maps on Fig [H We find that the best 
agreement is achieved when the symmetry of defects is 
described by characteristic tensor L with only non-zero 
elements L^x = Lyy = —2Lzz' The defects of this type 
tend to compress the surrounding lattice in x and y di- 
rections while elongating it in the remaining z direction 
or vice versa. The volume of unit cell tends to be pre- 
served. The elastic constants of morphotropic PZT single 
crystals, needed for calculations, have not yet been de- 
termined experimentally, but ceramics data (see Ref. l2l| 
and references therein) and theoretical estimations [22| 
are available. Pseudocubic elastic constants extracted 
from ceramics data [21] allow us to get a satisfactory 
qualitative description of DS distributions. However the 
best agreement we find with slightly changed constants 
cii = 135, Ci2 = 75 and C44 = 70 GPa. The correspond- 
ing comparison is presented in Fig. IH 

Despite apparent simplicity of our model we find it 
very reasonable. While we do not know the exact sym- 
metry of the host phase we may assume it to be close 
to cubic. First because the deviations from cubic struc- 
ture at least for compositions close to MPB are rather 
small. Secondly, after averaging over all possible orien- 
tations of monoclinic/rhombohedral domains the matrix 
will appear effectively cubic. The assumption of tetrago- 
nal symmetry of clusters is also well-founded since previ- 
ous studies indicate the signs of tetragonal phase in the 
MPB region. And, at last, this model perfectly repro- 
duces all the main peculiarities of the DS revealed by our 
experiments. It gives a zero-intensity transverse plane in 
high-symmetry (0 L) zones and a non-zero longitudinal 
scattering. Also it reproduces the observed q~'^ scatter- 
ing law. This combination could not be reproduced by 
purely form- factor based models such as the ones, pro- 
posed in Ref. HH. 
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FIG. 4: A comparison of the experimental 3-D distribution of diffuse scattering (left) and calculations based on the model of 
Huang scattering due to tetragonal defects (right). Pseudocubic elastic constants cn = 135, Ci2 = 75 and C44 = 70 GPa were 
used for calculations. 



In summary, in this Letter we reported complementary 
diffuse and inelastic X-ray scattering measurements on a 
morphotropic PZT single crystal that in tandem allow us 
to precisely point to the most probable microscopic or- 
ganization of the material below the morphotropic phase 
transition temperature Tmpb- When the temperature 
falls below Tmpb the tetragonal phase is not completely 
destroyed but remains in the form of local inclusions 
within the host phase of a different symmetry. We do not 
see any decline in the corresponding DS down to room 
temperature, which indicates a high stability of the pro- 
posed structurally heterogeneous state. 
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